On a new symmetry of the solutions of the wave 
equation in the background of a Kerr black hole 

Horst R. Beyer, Irina Craciun* 

Louisiana State University (LSU) 
Center for Computation and Technology (CCT), 
330 Johnston Hall 
& 

Department of Mathematics, 
148 LockettHall 

Baton Rouge, LA 70803, USA 1 
7th February 2008 



'Gratefully acknowledges financial support from CCT and the Department of Mathematics 
f E-mail: horst@cct.lsu.edu, Home Page: http://www.aei.mpg.de/_horst. 



1 



Abstract 



This short paper derives the constant of motion of a scalar field in the grav- 
itational field of a Kerr black hole which is associated to a Killing tensor of 
that space-time. In addition, there is found a related new symmetry operator 
S for the solutions of the wave equation in that background. That operator 
is a partial differential operator with a leading order time derivative of the 
first order that commutes with a normal form of the wave operator. That 
form is obtained by multiplication of the wave operator from the left with 
the reciprocal of the coefficient function of its second order time derivative. 
It is shown that S induces an operator S that commutes with the generator 
of time evolution in a formulation of the initial value problem for the wave 
equation in the setting of strongly continuous semigroups. 

1 Introduction 

In a remarkable paper, Carter discovered that the Hamilton-Jacobi equation is com- 
pletely separable for the Kerr space-time (M,g).[ll] This result was unexpected 
since (M, g) admits only a two-parameter group of isometries which by Noether's 
theorem leads to two independent conserved quantities for fields in a Kerr back- 
ground. 

The complete separability of the equation pointed to the existence of another less 
obvious 'symmetry' of (M,g). Indeed, this property could be traced back to the 
existence of a Killing tensor field K which is not induced by the isometries of 
(M, g).[48] It turned out that, together with the existence of two linearly inde- 
pendent Killing vector fields and the existence of a further Killing-Yano tensor 
field, this fact implies also the complete separability of the equations for Spin 
0, 1/2, 1 and 2 fields for Kerr space-time. [13, 31, 29, 30] On the other hand, all 
these equations have been separated before without reference to that tensor field 
which might be the reason that some consequences of the additional symmetry 
property of (M,g) are apparently not generally known.[ll, 14, 15, 16] 

Such a consequence is that to K there are also associated symmetry operators 1 
which commute with the differential operator denning those field equations. [12, 

1 Here and in the following, we call an operator a symmetry operator (for the solutions of) a field 
equation if it induces a map between classical solution spaces containing solutions with compo- 
nents belonging to a certain C*-class for some k e N*. 
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13, 39, 29, 30] For instance in the case of scalar fields, with respect to a local chart 
x = (x°, . . . , x 3 ) of M, such an operator is given by 1 

□ := -F=d a ^\K ab d h . (1.0.1) 
-\g\ 



It commutes with the wave operator 

1 

Rg\ 



i.e., 

(□□-□□)« = 

for all C 4 -functions u defined on the domain of x. In this, | | denotes the determi- 
nant function. 



As a consequence, it follows by Noether's theorem the existence of constants of 
motion that are associated to classical solutions of the field equations of sufficiently 
high order of differentiability. Considering the importance of constants of motion 
in the analysis of physical systems, these quantities could play an important role 
in the stability discussion of fields on a Ken - background where they have not been 
used so far.[l, 2, 3, 4, 9, 17, 19, 25, 23, 24, 27, 31, 34, 35, 37, 38, 42, 47, 46, 49, 52] 
Such potential importance is also suggested by the fact that in the description of 
classical point particles the constant of motion 

K(c,c) , (1.0.2) 

which is associated to K given by (2.0.4), is strictly positive on non-trivial light- 
like and time-like geodesies c were c denotes the tangent vector field along c. This 
is different from the constant of motion ('energy') corresponding to the ('Killing-') 
time coordinate vector field associated with Boyer-Lindquist coordinates that can 
assume negative values. The same is also true for scalar fields which considerably 
complicates their stability discussion in a Kerr background. It might be possible to 
simplify the last problem by taking into account the additional symmetry property 
of(M,g). 

1 In the following, Einstein's summation convention is adopted whereby a repeated index implies 
summation over all values of that index. In particular, Latin indices run from to 3. The symbol 
d a denotes the partial derivative operator in the x°-th coordinate direction for every a e {0, ... ,3}. 
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For these reasons, in this paper, we derive the constant of motion which is as- 
sociated to K for the case of the Klein-Gordon field in a Ken - background using 
Boyer-Lindquist coordinates. The derivation is an application of a more general 
discussion of Noether's theorem for PDEs that does not invoke a variational prin- 
ciple. That approach is essentially hidden in [12] and is briefly sketched in the 
Appendix. 

A drawback of the known symmetry operators that are associated to K is that 
their leading order time derivative is of the same or higher order than that of the 
field equations. [13, 31, 29, 30] This is different from the symmetry operators in- 
duced by Killing vector fields and presents an obstacle for their use in a rigorous 
functional analytic treatment of fields in a Kerr background. In such a setting 1 , 
the field equation is reformulated as a differential equation of the form 

u\t) = -Gu(t) 

for all tel. Here u is a map from the real numbers into a Banach ('state'-) space 
X, ' denotes the ordinary derivative of X-valued maps and G is a densely-defined, 
linear operator in X generating time evolution. See [4] for the details in the case of 
the Klein-Gordon equation in the background of a Kerr black hole. 

In such a context, it would be desirable that a symmetry operator would 'com- 
mute' in the following sense with G. It should lead on an intertwining operator 
2 S, i.e., a densely-defined linear operator that intertwines between the operators 
of the group of time evolution generated by A. [3 3, 6] Only a symmetry opera- 
tor whose leading order time derivative is of smaller order than that of the field 
equation leads in a simple way to a candidate for the restriction of an intertwining 
operator to a dense domain in X. 

In this paper, it will be shown that such a symmetry operator S indeed exists 
for the wave equation on a Kerr background in a normalized form. That form 
is obtained from the original geometric form by multiplication with the reciprocal 
of the coefficient of the leading order time derivative. Indeed, it is the last form 
which is the starting point for the reformulation of the wave equation in a func- 

1 For instance, see [6, 8, 18, 20, 21, 22, 26, 28, 32, 36, 43, 45, 50]. 

2 Such operators are also sometimes referred to as 'constraint' operators. 



4 



tional analytic setting. [4] As a consequence, differently to other known symmetry 
operators connected to K, that symmetry operator leads on a operator S which is 
defined on dense domain of X where it commutes with the infinitesimal generator 
of time evolution G given in [4]. Hence S satisfies the most basic requirement for 
an intertwining operator. It is 'time-independent' in the the sense that it does not 
contain time derivatives of the components of u. It is a candidate for the infinites- 
imal generator of a strongly continuous group or semigroup of symmetries on the 
data space X which are compatible with the time evolution. 



2 Constants of motion for the solutions of the Klein-Gordon 
equation in a Kerr background 

In Boyer-Lindquist coordinates 1 (t, r,0,<p) : O — > R 4 the Kerr metric g is given by 
[7] 

/ 2Mr\ 2Marsm 2 9 
g - 1 — \dt ® dt + (dt ® dtp + dip® dt) 



I . . / 2 2 2Ma 2 rsin 2 9\ , 



- — dr <8> dr - ~Ld& ® dQ - \ r ■ + a + j sin 6 dtp ® dtp 

where M is the mass, a e [0, M] is the rotational parameter, 

A := r 2 - 2Mr + a 2 , £ := r 2 + a 2 cos 2 6» , r + := M + Vm 2 - a 2 , 
Q := R x (r+, oo) x {-n, n) x (0, n) . 

As a reminder on the Ken - geometry, we give the following basic facts. The coor- 
dinate vector field d/dr becomes singular on the event horizon of Ken - space-time 
given by r = r+. The coordinate vector field d/dt is null on the ergosphere given 

by 

r - M + Vm 2 - a 2 cos 2 6 , 

it is space-like inside and time-like outside of ergosphere. Finally, the Kerr met- 
ric is globally hyperbolic outside the horizon and hence the Cauchy problem for 

1 In the following, the symbols t, r, 9, tp denote coordinate projections if not otherwise indicated. 
The domains of those projections will be obvious from the context. 
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the Klein-Gordon equation is well posed for data on any Cauchy surface. In the 
following, we consider only solutions of the Klein-Gordon equation of the form 



uit, r, 9, <p) - u(t, r, 9) e 



imsp 



for (t,r,6,<p) e £1 where m e Z. The reduced wave operator □„,, corresponding 
to □ by the replacement of d v through multiplication by the constant function of 
value im, is given by 



(r 2 + a 2 ) 2 2 • 2, 
■ a sin 6 



1 

m ■- ^ 

dgsinddg - m L \ '— - 

sin© 



d 2 + 



AimMar 



■ d t - d,-Ad r 



sin 2 © 



and the reduced Klein-Gordon equation for a field u of mass mo > by 

D m u + MqU = . 



(2.0.3) 



A Killing tensor K of the Kerr metric is given by 



K = a 2 1 



2Mrcos 2 © 



dt ® dt 



a sin 2 9 



[^(r 2 + a 2 ) + a 2 cos 2 © A] (dt ® dip + dtp® dt) 



- a 1 cos 2 © — dr®dr + r z I,d9®d9 
A 



sin 2 © 



+ 



pCr 2 + a 2 ) 2 + a 4 sin 2 © cos 2 © A] dip ® dip . (2.0.4) 



Hence it follows by (1.0.1) that the reduced operator □,„, corresponding to □ by 
the replacement of through multiplication by the constant function of value im, 
is given by 



r 2 sin 2 © + cos 2 ( 



(r 2 + a 2 ) 2 



d 2 



+2ima \r 2 + a 2 cos 2 © 



tt + a 



d t - a cos-9d r Ad r 



r 2 2 / r 2 a 4 cos 2 © 

+ -— do sin© da - m 



sin© 



sin 7 © 
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Note that this operator reduces in the special case of a Schwarzschild metric a = 
to 



Og S1IW Og + 



sin# 



sin 2 6> 



The eigenfunctions of the Sturm-Liouville operator that is induced by the last op- 
erator are used in the standard reduction of the wave equation in the gravitational 
field of a Schwarzschild black hole to a sequence of decoupled wave equations in 
one space and time dimension. 

A hyperbolic symmetry operator used later on is given by 



n - 2 _- .1 

□ m . — a. □ m □ m — 



2Ma 2 r(r 2 + a z )sin z 6» 2 



lima 99 99 99 r 2 + a 2 

(r 2 + a 2 ){A - a 2 sm 2 9)d t - a 2 sm 2 6d r Ad r : 

A smd 



dg sin^ dg 



+m 



, 1 r 2 + a 2 a 4 sin 2 ' 



sin 2 6> A 

In the following, □ m and ~u m are assumed to have C 2 (H,., C) as their domains where 

Q. r :- Ix (r + ,oo) x (0,tt) . 

Using the approach to Noether's theorem from Section 5.1 in the Appendix and 
after some calculation, it follows that the constant of motion corresponding to □ is 
given 



K{u) := -i J [j{u, □ u)] Q drdOdip 
Ja 

= 2n \ 



(2.0.5) 



/(r + ,oo)x(0,7T) 



sin# < 2ma 



+ .„ _ |2 
: \9qu\ 



+A\d r u\ + 



ir + a z 



\dgu\ 2 



m 2 I 



u\ 2 + -2— (A - a 2 sm 2 9)\u 2 \ 



+ 



Asin 2 6»'"' A 
+ia 2 sin 2 # A [(d r doii)* d r u — (d r u)*d r dou] 
(ji + a ~) 

+i [{dgd u}* dgu - {dgu)* dgdou] 
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+— 




2mla 2 Mr(r 2 + a 2 ) sin 2 6>] ((d u)*u - u*d u)} drdQ 



where for simplicity u e C 2 (O r ,C) is assumed as a solution of (2.0.3) such that 
u(t, •) has a compact support contained in 

Q„ := (r+,oo) x (0,;r) 

for every ? from some open interval of R. 

Since constants of motion can be combined into other such constants, we give 
also those corresponding to the Killing vector fields d t and d v , i.e., the energy E{u) 
and the z-component of the angular momentum L z (u). In this again, for simplic- 
ity, it is assumed that u e C 2 (Q r , C) is a solution of (2.0.3) such that u(t, •) has a 
compact support contained in Q. rs for every t from some open interval of R. 



3 A new symmetry of the solutions of the wave equation 
in the background of Kerr black hole 

For the derivation of such an operator, we first give a proof of the fact that □ m and 
□ m commute. This fact is a consequence of a more general result in [12]. For the 




and 




(r+,oo)x(0,n-) A 

-u*dou] + 4m 2 Mar\u\ 2 }drd6 . 



u 
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special case considered in this paper, a simpler direct proof is possible. Also, a 
similar proof will be given for the new symmetry operator. 



Theorem 3.1. Let u e C 4 (Q, r , C). Then 

LJ ffj LJ — I— I m I— I mil . 

Proof. Following [29], we define for every n eZ the linear operators D n , D+, L n , L+ 
C l (n r ,C) -> C(Q r ,C) by 



jd r + — [(r 2 + a 2 )d t + /ma J + 



2n(r - M) 



D> := jd r - ^ [(r 2 + a 2 )d f + /ma] + 2n(r ^ M) | 



L„v : 



m \ cos# 

3g + I —ia sin# <9 f H + 

sm6>/ sin# 

m \ cos# 



dg - I —ia sin# 3, + — — + n — 



sin#/ sin# 



for every v e C l (Q. r , C). In addition, we define D 2 , L 2 , L+ : C 2 (Q r , ( 
by 



C(Q r , C) 



D 2 v := (DjD+ +D|D )v 
2 2 

= — d r Ad r \(r 2 + a 2 ) 2 d 2 + limair 2 + a 2 )d t - m 2 a 2 ] v , 

A A 2 L J 

L 2 v := (LjL+ + L|L )v 



dg sin# dg + a sin~# df + limad, — — v 



sin 2 <9 



for every v e C 2 (Q r , C). As a consequence, it follows that 

□ m v - -— (AD 2 + L 2 ) v , v = -— (a 2 cos 2 AD 2 - r 2 L 2 ) • 
for all v e C 2 (n,., C) and, finally, that 

[□ m , □ ,„]« 



— a 2 cos 2 6>AD 2 , — L 2 
2E 2E 



w - 



— L z , — AD 

21 21 



9 



Aa 2 cos 2 6> 9 1 9 1 9 Aa 2 cos 2 6> , r 2 9 A 9 

D 2 —L 2 u L D m L — D 2 w 

2£ 2X TL TL TL TL 

A 9 r 2 9 A a 2 cos 2 # 9 A 9 r 2 9 

+ — D — L b - — D L u H D — L 2 w 

2E 2E 2Z 2X 21 2E 

1 9 Ao! 2 cos 2 6> 9 1 9 r 2 A 9 

L Dm L D u 

TL TL TL TL 

= — D 2 L 2 u - — L 2 AD 2 u = . 
41 4L 



Formally, the 'removal' of the highest order time derivative from □ m can be achieved 
as follows. For this, let u be a solution of 

n m u = 0. (3.0.6) 

The last equation can be solved for the second order time derivative of u. By 
substitution of that derivative into the expression u m u, there is obtained a partial 
operator S satisfying 

□ m u = S u 

with a leading order time derivative which is of the first order. That operator has 
an obvious extension to C 2 (Q. r , C) which we denote by the same symbol 5 . It is 
given by 

Su l - =- Pu , (3.0.7) 

(r 2 + a 2 )X + 2Ma 2 rsin 2 <9 



Pu := 



22 22 (r + a ) 
-2ima(r + a )Ld t - a sin 9 Ad r Ad r dg sin^e^ 

m 2 L „ 



+ — t-(Z + 2cT sin z 6») 
sin 

for every u e C 2 (O r , C). By construction, S is a symmetry operator for the so- 
lutions of (3.0.6) since it maps C 4 -solutions of (3.0.6) into C 2 -solutions. On the 
other hand, it is not difficult to see that S does not commute with □,„. The com- 
mutator is given in formula 6.0.4 in the Appendix. Somewhat surprisingly, it turns 
out that S commutes with the operator 

AL 9 MmMar 

dt + 9- d t 



(r 2 + a 2 )L + 2Ma 2 r sin 2 6> (r 2 + a 2 )L + 2Ma 2 r sin 2 6* 
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1 



(r 2 + a 2 )T, + 2Ma 2 rsm 2 9 



sin# 



+ -^-(A-a 2 sin 2 (9) 
sin 2 <9 



Theorem 3.2. Let u € C A {Q. r , C). Then 
AE 



□ m S U - S 



AE 



(r 2 + a 2 )E + 2Ma 2 r sin 2 6» (r 2 + a 2 )!, + 2Ma 2 r sin 2 6» 

Proof. To simplify notation in the following, we define 

(r 2 + a 2 )Z + 2Ma 2 r sin 2 6> (r 2 + a 2 ) 2 , . , 
I := - ■ a sin 8 . 



Then it follows from (3.0.6) that 



AS 



(r 2 + a 2 )E + 2Ma 2 rsin 2 6> 



□mV 



= _±(ad 2 + l 2 V 

21 V ; 



for every v e C 2 (Q r , C). Further, by a straightforward calculation, it follows that 
p v = -- [a 2 sin 2 A 2 D 2 + (r 2 + a 2 ) 2 L 2 ] v 



and hence that 



1 

Sv = -- 



1 



2 (r 2 +a 2 )E + 2M£ Z 2 rsin 2 <9 



[a 2 sin 2 A 2 D 2 + (r 2 + a 2 ) 2 L 2 ] v 



1 

21 



( r 2 _j_ fl 2\2 

a 2 sin 2 AD 2 + L 2 

A 



for every v e C 2 (fi,-, C). Finally, it follows that 



[5, EE 1 □„,]« = -{- 



a 2 sin 2 0AD 2 + (r2+a2)2 L 2 



1 

4E 



-(AD 2 + L 2 ) m 

a 2 sin 2 AD 2 = L 2 - AD 2 = { — + ^ L 2 
E E A 
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(r z + a 2 ) 2 9 A 7 9I9 ? 9 
L 2 -D 2 - L — a 2 sin 6 AD 

4SL A I S 

A „ 2 a 2 sin 2 6> - ((r 2 + a 2 ) 2 /A) 2 



= — rr = '-^ u-u 

A 2 ((r 2 + a 2 ) 2 1 A) - a 2 sin 2 (9 2 
41 I 

= -— D 2 L 2 m + —= L 2 D 2 u — . 
4X 4X 



□ 



Note that the constants of motion associated to a solution u of the reduced wave 
equation (3.0.6) and the operators n m ,a m , S are, apart from a sign, the same since 

In [4], it was proved the well-posedness of the initial value problem for the re- 
duced Klein-Gordon in a Kerr background in the setting of strongly continuous 
semigroups of operators, the domain of the corresponding infinitesimal generator 
G was investigated and its resolvent given in terms of spheroidal functions. Fur- 
ther, the stability of the solutions was proved if the mass of the field exceeds a 
certain combination of the 'magnetic quantum number' m e Z relative to the rota- 
tion axis, the mass and the rotational parameter of the black hole. For details refer 
to [4]. In particular for (it, v) e C 2 (Q„, C) x C (Q„. Q 

G(u, v) - (-v, (A + C)u + iBv) 



where 



(A + C)u = 



1 



(r 2 +c? 2 )X + 2Ma 2 rsin 2 6» 



-Ad r Ad r 



sin# 



■ da sin# da 



+—, T-(A-a z sin z 6>) 



sin 2 <9 



Bv = 



AmMar 



v . 



(r 2 +a 2 )T, + 2Ma 2 rsm 2 G 

It is not difficult to see that in the setting of [4], S induces the operator S defined 
by 

* 1 9 9 9 9 

S(u, v) = — = (P s u — 2ima(r + a )Ev, limair + cr)£(A + C)u 
AS 
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+ P s v - 2ma(r 2 + a 2 )Sfiv) 



for all u, v € Cq(Q,-. s , C) where 



P s u := 



a 2 sin 2 9 Ad,-Ad r 



(r 2 + a 2 ) 2 



sin9 



do sin# dg 



m 2 S 9 9 

+ — — (S + 2ct 2 sin 2 6») 
sin" 9 



it . 



Then the following holds 

Theorem 3.3. Let (a, v) e C*(Q. rs , C) x C 2 (Q„, C). Then 

SG(u, v) = GS(u, v) . 
Proof. First, it follows from Theorem 3.2 that 



-lima 



(r 2 + a 2 )S 



—PA + C 
AS 



AS 

W2 = 



+ 1 



1 

—P,B 
AS 



w\ = , 



for all wi € C 2 (Q„, C) and w 2 e Cj(n„, C). Further, it follows 
SG(u, v) 

- -L(-2zma(r 2 + a 2 )S(A + C)u - Pv + 2ma(r 2 + a 2 )I,Bv, 
AS 

P(A + On - 2ma(r 2 + a 2 )SB(A + C)w 

- 2/ma(r 2 + a 2 )S(A + C)v + iPBv - 2ima(r 2 + a 2 )LB 2 v) 
GS (u,v) 

= X=.(-2ima(r 2 + a 2 )S(A + C)u - Pv + 2ma(r 2 + a 2 )SBv, 
AS 

AS(A + Q—Pu - 2ma{r 2 + a 2 )l,B(A + C)u 
AS 

- 2/maAS(A + C)- ^—v + iBPv 

AS 

- 2ima(r 2 + a 2 )SB 2 v) 

and hence by (3.0.10) the validity of (3.0.9). 
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4 Discussion and open problems 



The results open up a number of interesting new questions for future research. 

It would be interesting to know whether and, if applicable, when the new constant 
of motion K{u) 2.0.5 (or some combination with the other constants of motion) 
has a definite sign like its classical counterpart (1 .0.2). For instance, a definite sign 
of K(u) in the case of a vanishing mass of the field would open up the road to a 
simpler and independent proof of the stability of the solutions of the wave equa- 
tion in a KeiT background. [23] In addition, it could lead to a sharper estimate on 
a possible onset of instability in the case of a non- vanishing mass of the field. [4, 35] 

Connected to the new symmetry operator S (3.0.7) are technical questions whether 
the closure of the induced operator S is an intertwining operator and/or is the in- 
finitesimal generator of a strongly continuous group or semigroup. Also, with the 
help of the new symmetry operator S , it should be possible to further reduce the 
wave equation (3.0.6), in a proper sense, to an equation in 'one space dimension'. 
For this, a deeper study of the properties of S is necessary. In addition, it would be 
interesting to explore whether there can be made a connection of the results here 
to those of [49]. Finally, it remains the question whether such a symmetry operator 
can also be found for a non-vanishing mass of the scalar field. 

It is to expect that the results of the paper can be generalized to the Dirac, Maxwell 
and the linearized Einstein's field equation around the Kerr metric. If that turns out 
be the case, it would be interesting to investigate their potential for finding a gen- 
eralization to fields on a Kerr background of the Regge-Wheeler-Zerilli-Moncrief 
decomposition of fields on a Schwarzschild background. [44, 51, 40] 

5 Appendix 

5.1 An approach to Noether's theorem not invoking a variational prin- 
ciple 

In this Section, we briefly give an approach to the derivation of conservation laws 
from transformations mapping classical solution spaces of a linear scalar PDE into 
one another. The approach is essentially hidden in [12] and is not invoking a vari- 
ational principle. The derivation of the constants of motion in Section 2 is an 
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application of this approach. 



For this, we consider the equation 

1 _d_ 

P dx a 



Pg 



du 
dx h 



+ Vu = 



(5.1.1) 



where n e N*, f2 is a non-empty open subset of R" +1 , g" h e C^QjR) for all 
a,be{0,..., n] such that g ah = g ba for all a, b € {0, . . . , n], p e C 1 ^, R) such that 
Ranp c (0, oo), V e C(Q, R) and u e C 2 (Q, C). 



For solutions u, v of (5.1.1), it can be derived a corresponding conservation law. 
First, it follows that 



P dx a 



dx b 



dxbj 



where * denotes complex conjugation. The difference of these equations gives 



P OX a 



PS 



ab 



dv \ v d 
dx b J p dx a 



ab 



dx a 



ab 



dv 



1 d 

P dx a 



U pg 



„ab I ..* 



d 



dx b J dx a 



pg 



vpg 



du* 
dx h 

ab^£_ 

dx b 



dv 



and hence 



Pg « 



- v 



du* \ 



dxb dxb J 



-d a p[j(u,v)] a = 
P 



where the vector field j(u, v) e C 1 (Q., R" +1 ) is defined by 



dv 



du" 



[;( „,v)r= s -|„---v- 



for all a e {0, ... , n}. If F : C k (Q, C) C 2 (0, C) for some it 6 {2, 3, ... } maps 
C^-solutions of (5.1.1) into C 2 -solutions of (5.1.1), then it follows the constraint 

-d a p[j(u,F(u))] a = 
P 
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for every C^-solution u of (5.1.1). 



As a short example, we consider the special case of a wave equation in flat space- 
time 

3 2 u 



dt 2 



where V e C(Q, R) is such that 



- Aw + Vu = 



^=0 

dt 



and u is assumed as real-valued. If u is a C -solution of (5.1.2), then 



(5.1.2) 



a t&u 

~ diVdfl 



d 2 du . du 



du 



Au + Vu\ = — A — + V — 

dt 2 dt dt dt 



and hence du/dt is a C 2 -solution of (5.1.2). As a consequence, 



d_ 

at 



du 
1 1 u, — 

dt 



-y— 



du 
j\u, — 

of 



Further, 



where 



du 

J \ U >-dt 



— uAu 



d 2 u ldu\ (du 



du , 

a 1 



= V (uVu) 



— I +\Vu\ z + Vu l 
dt 



V (wVw) - 2 e(u) 



eiu) := - 



^) + |Vw| 2 + V M 2 
of / 



and V denotes the gradient in the spatial coordinates. In addition, 



du 



— —\u 



d 2 u du du \ 



dxkdt dt dxk] 



Hence it follows that 



d I du du 

— tV(«V«)-2e(i«)]=Vl«V— -— V»| . 
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The last implies that 



!<- 2£ <"»=- 2v w va 



and, finally, the conservation law 

de{u) 
dt 



v,'|v„ 



(5.1.3) 



As is well-known, in physical applications, e{u) is called the energy density and 

du _ 

is the energy flux density associated with u. We note that (5.1.3) is true also for 
C 2 -solutions of (5.1.2) since for such an u 



dt \dt dt2 



-) + \Vu\ 2 + Vu 2 
dt J 



V — Vm 

dt 



du d 2 u d w du (_ d \ _ s du k 

= + (Vm) • V— u + Vu V— u ■ (Vm) Au 

dt dt 2 dt dt \ dt j dt 

du t d 2 u 

= *U- A " +l "' l=0 ' 



6 The commutator of the wave operator and S 

In the following we give without proof the commutator of the wave operator and 
the new symmetry operator S : 



[n,S]u = (Ci + C 2 d r + C 3 d e ) □ u 



(6.0.4) 



where 



T 2 



2Ma' 



E(AX) 



-{sin 2 0(rir 2 -A£Ar 3 ) 



+2r(r 2 + a 2 ) 3 [a£ (2cos 2 6» - sin 2 6») + 4a 2 Asin 2 0cos 2 fl]) , 
:= (r 4 -a 4 )I + 2a 2 r 2 Asin 2 6» , 
:= 2 [(r - M)A£ + 2r(r 2 + a 2 ) A - 2M(r 4 - a 4 )] , 
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r 3 := 2 [r(r 4 - a 4 ) + 2r 3 Z + 2a 2 rA sin 2 6> + 2a 2 r 2 (r - M) sin 2 6»] , 
4Ma 2 A sin 2 6» [(r 4 - a 4 )X + 2a 2 r 2 A sin 2 #] 



C 2 := 
C 3 := 



i(A£) 2 

%Ma 2 r(r 2 + a 2 ) 3 sin cos 
S(AX) 2 



and w € C 4 (Q r , C). Note that the form of the commutator (6.0.4) generalizes that 
in the definition of a symmetry operator given in [39] to symmetry operators of the 
second order. 
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